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Abstract — This paper considers Gaussian relay networks wliere 
a source transmits a message to a sink terminal with the help of 
one or more relay nodes. The relays work in half-duplex mode, in 
the sense that they can not transmit and receive at the same time. 
For the case of one relay, the generalized Degrees-of-Freedom 
is characterized first and then it is shown that capacity can 
be achieved to within a constant gap regardless of the actual 
value of the channel parameters. Different achievable schemes 
are presented with either deterministic or random switch for the 
relay node. It is shown that random switch in general achieves 
higher rates than deterministic switch. For the case of K relays, 
it is shown that the generalized Degrees-of-Freedom can be 
obtained by solving a linear program and that capacity can be 
achieved to within a constant gap of A721og(4A"). This gap may 
be further decreased by considering more structured networks 
such as, for example, the diamond network. 

Index Terms — Relay Channel, Generalized Degrees-of- 
Freedom, Capacity to within a Constant Gap, Inner bound. 
Outer bound, Half-duplex. 



I. Introduction 

The performance of wireless systems can be enhanced by 
enabling cooperation between the wireless nodes. The simplest 
form of cooperation is modeled by the Relay Channel (RC) 
where a source terminal communicates to a destination with 
the help of a relay node. In this multi-hop system the relay 
helps to increase the coverage and the throughput of the 
network. 

Relays employed in practical wireless networks can be 
classified into two categories: Full-Duplex (FD) and Half- 
Duplex (HD). In the former case the relays transmit and 
receive simultaneously; in the latter case the relays can either 
transmit or receive at any given time, but not both. There are 
some relatively expensive relay devices which work in FD 
mode, normally used in military communications. However 
FD relaying in wireless networks has practical restrictions such 
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as self-interference, which make the implementation of decod- 
ing algorithm challenging. As a result HD relaying proves to 
be a more practical technology with its relatively simple signal 
processing. Thus it is more realistic to assume that the relay 
operates in HD mode either in Frequency Division Duplexing 
(FDD) or Time Division Duplexing (TDD). In FDD, the 
relay uses one frequency band to transmit and another one 
to receive; in TDD, the relay listens for a fraction 7 G [0, 1] 
of time and then transmits in the remaining 1 — 7 time fraction. 
From an application point of view, the HD model fits future 
4G network with relays [IJ, where the relay communicates 
over-the-air with the source, which is called the Donor-eNB. 
We keep our focus on deployment scenarios where the relay 
works in TDD HD mode. 

In this work we concentrate on the HD relay networks, 
where the relays transmit and receive in different time slots. 
HD relaying has received considerable attention lately, as 
summarized next. 

A. Related Work 

a) Single Relay Networks: The RC was first introduced 
by van der Meulen [2] and then thoroughly studied by Cover 
and El Gamal Q. In (jSj the authors study the general memo- 
ryless RC, derive inner and outer bounds on the capacity and 
establish the capacity for some classes of RCs. The proposed 
outer bound is now known as the max-flow min-cut outer 
bound, or cut-set for short, which can be extended to more 
general memoryless networks |4|. Two achievable relaying 
strategies were proposed in |3|, whose combination is still 
the largest known achievable rate for a general RC, namely 
Decode-and-Forward (DF) and Compress-and-Forward (CF). 
In DF, the relay fully decodes the message sent by the 
source and then coherently cooperates with the source to 
communicate this information to the destination. In CF, the 
relay does not attempt to recover the source message, but it 
just compresses the information received and then sends it to 
the destination. The capacity of the general memoryless RC 
is known for some special classes, in particular for degraded 
RCs, reversely degraded RCs and semi-deterministic RC |3J. 

The HD-RC was studied by Host-Madsen in 0. Here 
the author derives both an upper and a lower bound on the 
capacity. The former is based on the cut-set arguments, the 
latter exploits the Partial-Decode-and-Forward (PDF) strategy 
where the relay only decodes part of the message sent by the 
source. Host-Madsen considers the transmit/listen state of the 
relay as fixed and therefore known a priori to all nodes. 

In [6 1, Kramer shows that larger rates can be achieved by 
using a random transmit/listen switch strategy at the relay. In 
this way, the source and the relay can harness the randomness 
that lies in the switch in order to transmit extra information. An 
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important observation of ||6l is that there is no need to develop 
a separate theory for memoryless networks with HD nodes as 
the HD constraints can be incorporated into the memoryless 
FD framework. In this work we shall adopt this approach in 
deriving outer and inner bounds for HD relay networks. 

b) Multiple Relay Networks: The pioneering work of 
Cover and El Gamal |3J has been extended to networks with 
multiple relays. In Q Kramer et al. proposed several inner 
and outer bounds as a generalization of DF, CF and the cut- 
set bound. It was shown that DF achieves the ergodic capacity 
of a wireless Gaussian network with phase fading if phase 
information is available only locally and the relays are close 
to the source node. 

The exact characterization of the capacity region of a 
general memoryless network is challenging. Recently it has 
been advocated that progress can be made towards under- 
standing the capacity by showing that achievable strategies 
are provably "close" to (easily computable) outer bounds H). 
As an example, Etkin et al. characterized the capacity of the 
Gaussian Interference Channel to within 1 bit regardless of the 
system parameters 191. In [lO], the authors study unicast and 
multicast Gaussian relay networks with K nodes and show that 
capacity can be achieved to within X^a^i 5 min{Mfc, A^^;} bits 
with quantize-remap-and-fonvard (QMF), where and Nf^ 
are the number of transmit and receive antennas, respectively, 
of node fc € [1 : X]. Interestingly, the result is valid for static 
and ergodic fading networks where the nodes operate either 
in FD mode or in HD mode with deterministic listen/transmit 
schedule for the relays. For single antenna systems, Lim et al. 
in 1 11 1 recently showed that this bK bits gap can be reduced 
to Q.QK bits for FD relay networks with noisy network coding 
(NNC). Both QMF and NNC ai-e network extensions of CF. 

The gap characterization of UJJ is valid for any multi relay 
network but linear in the number of nodes in the network, 
which could be a too coarse capacity characterization for 
networks with a large number of nodes. Tighter gaps can 
be obtained for more structured networks. For example, the 
diamond network model was first proposed in fT2l. A diamond 
network consists of a source, a destination and K — 2 relays. 
The source and the destination cannot communicate directly 
and the relays cannot communicate among themselves. In 
other words, a general Gaussian multi relay network with K 
nodes is characterized by K{K — 1) channel gains, while a 
diamond network only has 2{K — 2) non zero channel gains. 
In [12] the case of two relays was studied for which an achiev- 
able region based on time sharing between DF and amplify- 
and-forward (AF) was proposed. The capacity of a general FD 
diamond network is known to within 2\og{K — 1) bits ifTSl . 
If in addition the FD diamond network is symmetric, that is, 
all source-relay links are equal and all relay-destination links 
are equal, the gap is less than 2 bits for any K lfl4ll . HD 
diamond networks have been studied as well, albeit only for 
deterministic switch for the relays. In a HD diamond network 
with {K — 2) relays, there are 2^^^ possible combinations 
of listening and transmit states, since each relay can either 
transmit or receive. For the case of {K — 2) = 2 relays, ifTSll 
shows that out of 2^^^ = 4 possible states only {K — 1) = 3 
states suffice to achieve capacity to within less than 4 bits. 



Their achievable scheme is a clever extension of the two-hop 
DF strategy of [16|. It is interesting to note that fisl derived 
closed-form expressions for the fractions of time the relays are 
active and tight outer bounds based on the dual of the linear 
program (LP) associated with the classical cut-set bound. 
Extensions of these ideas to more than two relays appear 
difficult due to the combinatorial structure of the problem. 

Inspired by flSl, the authors in ||T71 showed that for a very 
specific HD diamond network with {K — 2) = 3 relays, {K — 
1) = 4 states out of 2^^^^ = 8 are active in the cut-set outer 
bound. In the same work, it was verified numerically that for a 
general HD diamond network with {K—2) < 7 relays, (if — 1) 
states suffice for the cut-set upper bound and it is conjectured 
that the same holds for any number of relays. We remark that 
in 1 17 1 only the cut-set upper bound was considered; moreover 
only the case of deterministic switch and per-symbol power 
constraint was considered. 

Multi relay networks were also studied in [18 ,1 where the 
authors determine numerically the optimal fractions of time 
each relay transmits/receives with DF through an iterative 
algorithm. Also in this case the relays use deterministic switch. 

B. Contributions 

In this work we focus on the HD relay networks. The exact 
capacity of this channel is unknown. In this paper we make 
progress toward determining its capacity by giving a constant 
gap result for any Gaussian network with random switch. Our 
main contribution can be summarized as follows: 

1) We determine the generalized Degrees-of-Freedom 
(gDoF) of the HD RC with a single relay. We identify 
three schemes that achieve the gDoF upper bound. The 
simplest one is inspired by the Linear Deterministic 
Approximation (LDA) of the Gaussian noise channel 
at high SNR |10|; it uses superposition coding at the 
source, DF at the relay and stripping decoding at both 
the relay and destination; we note that neither power al- 
location nor backward decoding is required at the nodes. 
The second and third schemes use more sophisticated 
coding techniques and are based on PDF and NNC 
strategies |4|. 

2) We prove that the three schemes above achieve the 
capacity to within a constant gap regardless of the 
channel parameters. We consider both deterministic and 
random switch for the relay. Thus in the second case the 
relay harnesses the randomness that lies in the switch to 
achieve larger rates and therefore smaller gaps from the 
cut-set upper bound. 

3) We prove that PDF with random switch is optimal for 
a Gaussian diamond network with one relay, i.e., it 
achieves the capacity, even though we were not able 
to determine the capacity achieving input distribution. 

4) We determine the capacity of the noiseless LDA channel. 
In particular we show that random switch and non- 
uniform inputs at the relay are optimal. 

5) For HD networks with K nodes (HD-MRC), all 
of which employ random switch, we prove that 
NNC achieves the cut-set outer bound to within 
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Kjl log(4iir) bits. For diamond relay networks assum- 
ing that the conjecture in flTl holds for any K, the gap 
can be reduced to 51og(A'). 

C. Paper Organization 

The rest of the paper is organized as follows. Section 
im describes the channel model and states our main result. 



Section III derives the gDoF upper-bound based on the cut- 
set bound. Section |IV] provides a lower bound based on the 
PDF strategy. Section|V]highlights a motivating example based 
on the LDA and provides a simple achievable based on this 



approach. Sections VI and VII are devoted to the analytical and 



numerical proofs that the capacity for a single relay network is 



achievable to within a constant gap, respectively. Section VIII 
considers networks with multiple relays and it shows that the 
gDoF can be computed by solving a linear program and that 
NNC achieves the capacity to within a constant gap. Section 
IX concludes the paper. 



II. Single Relay Networks: System Model 
A. General Memoryless Relay Channel 

A RC consists of two input alphabets {Xs,Xr), two output 
alphabets (3^r,3^d) and a transition probability PY.r.Ya\x,.Xr- 
The source sends symbols from Xs, the relay receives symbols 
in and sends symbols from AV, and the destination receives 
symbols in The source has a message £ [1 : 2^^] for 
the destination where N denotes the codeword length and R 
the transmission rate in bits per channel use At time i, 
i £ \1 : N], the source maps its message W into a channel 
input symbol Xg i{W) and the relay maps its past channel 
observations into a channel input symbol Xr/i{Y^^^). The 
channel is assumed to be memoryless, that is, the following 
Markov chain holds for all i e [1 : iV] 

At time N, the destination makes an estimate VF(y"j^) of the 
message W based on all its channel observations Y^ . A rate 
R is said to be e-achievable if V\W ^ W] < e for some 
e G [0, 1]. The capacity is the largest nonnegative rate that is 
e-achievable for any e > 0. 

We note that half-duplex channels are a special case of the 
memoryless full-duplex framework in the following sense IS): 
let the channel input of the relay be the pair {Xr, Sr), where 
Xr e Xr as before and Sr G {0, 1} is the state random 
variable that indicates whether the relay is in receive-mode 
(Sr ~ 0) or in transmit-mode (Sr = 1). The memoryless 
channel transition probability is defined as 



P. 



Yr^,Ya\X^,Xr,Sr=0 



R 



(0) 



Py 



Y^,Ya\Xs,Sr=0 

,(1) 



.R 



(1) 



Yr,Ya\X,,X^,Sr = l - - Ya\X,,X^,Sr = l^ y,|S, = l' 

that is, when the relay is in receive-mode (Sr = 0) the 
outputs Yr, Yd are independent of Xr and when the relay is in 
transmit-mode (Sr = 1) the relay output Yr is independent of 
everything else. In other words, the (still memoryless) channel 
is now specified by the two transition probabilities one for each 
mode of operation fS). 



B. The Gaussian Half-Duplex RC 

We consider a single-antenna complex-valued Gaussian 
Half-Duplex Relay Channel (G-HD-RC), shown in Fig. [T] 
where the inputs are subject to an average power constraint, 
described by the input/output relationship 

Yr = ^/CX, (1 - Sr) + Zr, (la) 

Yd = \fSXs + e^^VlXr Sr + Zd, Oe M, (lb) 

where the channel gains C, S, I are constant and therefore 
known to all terminals. Without loss of generality we can 
assume that the relay output Yr does not contain the relay 
input Xr because the relay node can always subtract Xr 
from Yr. Moreover, since a node can compensate for the 
phase of one of its channel gains, we can assume without 
loss of generality that the channel gains from the source to 
the other two terminals are real-valued and nonnegative. The 
channel inputs are subject to unitary average power constraints 
without loss of generality, i.e., E[|X„p] < 1, u G {s,r}. The 
'switch' random variable Sr is binary. The noises Zd, Zr are 
assumed to be zero-mean jointly Gaussian and with unit power 
without loss of generality. In particular (but not without loss 
of generality) in this work we assume that Zd is independent 
of Zr. In the following we will only consider G-HD-RC for 
which C > and / > in (|2]i, since for either C = or 
/ = the relay is disconnected from either the source or the 
destination, respectively, so the channel reduces to a point-to- 
point channel with capacity log(l + S). 

The capacity of the channel in ([T]) is unknown. Here we 
make progress toward determining its capacity by establishing 
its gDoF, i.e., an exact capacity characterization in the limit 
for infinite SNR |9|, and its capacity to within a constant gap 
at any finite SNR. Consider SNR > and the parameterization 



S 
I 

c 



= SNR^'^'', source-destination link, (2a) 

= S N R^'^'' , relay-destination link, (2b) 

= SNR^'^'. source-relay link, (2c) 
1. We define: 



(^(HD-RC) 



for some (/3sd, Ad, Ar) e 
Definition 1. The gDoF is 

j(HD-RC) ijj^ 

SNR^+oo log(l + SNR) ' 

where cCHD-RC) ^/^^ capacity of the G-HD-RC. 

Definition 2. The capacity (7^^^^^^) is said to be known to 
within b bits if one can show rates i?^™^ and such that 

^(in) ^ (^(HD-RC) < ^(out) < ^(in) + blog(2). 

Our main result for single relay networks can be summa- 
rized as 

Theorem 1. The gDoF of the G-HD-RC is given by ([3| at the 
top of next page and the cut-set upper bound is achieved to 
within the following number of bits 



Logarithms are in base 2. 



Achievable scheme 


LDA 


NNC 


PDF 


analytical gap 
numerical gap 
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1.59 


1.61 
1.52 


1 
1 
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^(HD-RC) 




for /3sr > /3sd, 
otherwise. 



/3rd > A 



sd 



(3) 



where LDA is a very simple achievable scheme inspired by the 
linear deterministic approximation of the G-HD-RC at high 
SNR, PDF is partial-decode-and-forward and NNC is noisy- 
network-coding, or compress-and-forward. 

Sections [lIll|VII| are devoted to the proof of Theorem [T| 

Remark 1. The gDoF of the Gaussian Full-Duplex Relay 
Channel (G-FD-RC) is 

j(FD-RC) ^ ^ min{[/3,, - Ad] + , [13rd ~ /3sd] + }, (4) 

and its capacity (7^^^^^^) is achievable to within 1 bit 
per channel use [10]. We notice that HD achieves the same 
gDoF of FD if mhi{/3j-d, Pst} < /3sd, in which case the RC 
behaves gDoF-wise like a point-to-point channel from TX to 
RX with gDoF given by /Sgd- In either FD or HD the gDoF 
has a 'routing ' interpretation UOV : if the weakest link from 
the source to the destination through the relay is smaller 
than the direct link from the source to the destination then 
direct transmission is optimal and the relay can be kept silent, 
otherwise it is optimal to communicate with the help of the 
relay. 



III. Single Relay Networks: Upper Bound 

This section is devoted to the proof of a number of upper 
bounds that we shall use for the converse part of Theorem [T] 
From the cut-set bound we have: 

Proposition 1. The capacity of the G-HD-RC is upper 
bounded as in (|5]l, (|6]l and (|7]l at the top of next page where 

• in (|5]l." the distribution g is the one that maxi- 
mizes the cut-set upper bound, 

• in (|6]l." the parameter 7 := V[Sr = 0] G [0, 1] represents 
the fraction of time the relay node listens, Hi"/) is the 
binary entropy function defined as 

m-f) :=-7log(7)-(l-7)log(l-7), (12) 

the maximization is over the set 



7e [0,1], 
l«i| < 1, 

(-Ps.O, -PsJ, -Pr.O, -Fr,l) G 



(13) 
(14) 



: iPu^o + (1 - l)PuS < 1, u G {s, r}, (15) 

and the mutual informations Ii, . . . , are defined as 

/i :=log(l + 5 P,,o), (16) 

h ■■= log (l + SP,,i+/P,,i + 2|ai|v/SP,,i/P,,i) , (17) 
h:=log{l + {C + S)Psfi), (18) 
/4:=log(l + (l-|ain^P,^i), (19) 



in the terms hi and 62 <2re defined as 

log(l + (V7 + \/5)2) 

61 := > 1 since I 

log (1 + S) 

log (1 + C + 5) 



log (1 + 5) 



> 0, (20) 
> 1 since C > 0. (21) 



Proof The proof and the definitions of the above quan- 
tities can be found in Appendix |A] ■ 
The upper bound in (|5| will be used to prove that PDF with 
random switch achieves capacity to within 1 bit, the one in (|6]) 
to prove that PDF with deterministic switch also achieves 
capacity to within 1 bit and for numerical evaluations (since 
we do not know the distribution P^^ that maximizes the 
cut-set upper bound in Q), and the one in (|7]) for analytical 
computations such as the derivation of the gDoF 

Proposition 2. The gDoF of the G-HD-RC is upper bounded 
by the right hand side of (|3]l. 

Proof: The proof can be found in Appendix IB] ■ 



IV. Single Relay Networks: Lower Bounds based 
ON partial-decode-and-forward 

This section is devoted to the proof of a number of lower 
bounds that we shall use for the direct part of Theorem [T] 
From the achievable rate with PDF we have: 

Proposition 3. The capacity of the G-HD-RC is lower 
bounded as in (|8]l, (|9]) and ( |10| l at the top of next page where 

• in we fix the input Pu.Xs.Xr,Sr evaluate the PDF 
lower bound; in particular we set Pxs,Xr,Sr- 

same distribution that maximizes the cut-set upper bound 
in (|5]l and we choose either U — or U = X^Sr + 

Xs{l-Sr). 

• in (|9]l." the parameter 7 :— V[Sr = 0] G [0, 1] represents 
the fraction of time the relay node listens, the maximiza- 
tion is over the set ( |13| l-( |15| l as for the cut-set upper 
bound in (|6]), the mutual informations /a , . . . , /§ are 

h := h in ([16]), (22) 

h ■■= h in ([nil, (23) 

h log (1 + max{C, 5}P«,o) < ^3 in ([18), (24) 

4 h in ([g, (25) 



and I, 



(PDF) _ 





I{Sr]Y^) is computed from the density 



(26) 



with t G C,vq = exp(/5), vi = exp(/g). 
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^(HD-RC) < j^i^^ I^I{X,,Xr,Sr;Yd),IiXs;Yr,Yd\Xr,Sr)} 

lin |-H(7) + 7/1 + (1 - 7)/2, 7^3 + (1 - 7)-f4} 



< max min ■ 



r 



(CS-HD) 



< 21og(2) + log(l + 5) 1 + 



(61 - l)(b2 - 1) 
ih - 1) + {b2 - 1) 



^(HD-RC) > j^^ij^ I y^)^ 

/(C/; Sr) + I{X,;Yd\Xr, Sr,U)} 

^(HD-RC) > j^ij^ |^(PDF) ^ ^ _ ^ _ ^^^^1 ^(PDF 

(Ci - 1)(C2 - 1) 
(Cl - 1) + (C2 - 1) 



HD) 



> log(l + S') 1 + 



(5) 
(6) 
(7) 

(8) 
(9) 
(10) 



log 1+ ITS 



log (1 + 1^5) + flog -log (l + T 



log 1 



C 



1 + 5 



-log 1 



S 



1 + 5 



T + 



(11) 



in P0| l; f/ie ferms Ci ant/ C2 are 

log (1 + / + g) 

log (1 + 5) 
log (1 + niax{C, 5}) 



Cl 



> 1 since / > 0, 



(27) 



C2 



> 1 since C > 0. (28) 



log (1 + 5) 

Proof: The proof can be found in Appendix |C] ■ 
The lower bound in ([8]) will be compared to the upper 
bound in (|5]l to prove that PDF with random switch achieves 
capacity to within 1 bit, the one in (|9]l with the one in (|6]l 
to prove that PDF with deterministic switch also achieves 
capacity to within 1 bit and for numerical evaluations, and 
the one in ( [TO] i for analytical computations such as evaluation 
of the achievable gDoF. 

Proposition 4. The gDoF of the G-HD-RC is lower bounded 
by the right hand side of (|3]l. 

Proof: The proof can be found in Appendix [D] ■ 
Propositions |2] and |4] show that the gDoF for the G-HD-RC 
is given by ([3]l. 

V. Single Relay Networks: A Simple Achievable 
Strategy 

In this section we propose a very simple achievable scheme 
that is gDoF optimal, that achieves capacity to within 3 bits 
and that can be implemented in practical HD relay networks. 
In Section IV-AI we describe a deterministic-switch achievable 
strategy for the Linear Deterministic Approximation (LDA) of 



the G-HD-RC at high SNR which we mimic in Section |V-B 
to derive an achievable rate for the G-HD-RC at any SNR. 
This achievable scheme is referred to as the LDA-strategy, or 
LDA for short. The main result of this section is: 

Proposition 5. The capacity of the G-HD-RC is lower 
bounded as in \\ \\ at the top of this page. 



The rest of the section is devoted to the proof of Proposi- 
tion |5] Before we provide the details of the scheme, we point 
out three important practical aspects of this scheme that are 
worth noticing: 

1) the destination does not use backward decoding, which 
simplifies the decoding procedure and incurs no delay, 

2) the destination uses successive decoding, which is sim- 
pler than joint decoding, and 

3) no power allocation is applied at the source or at the 
relay, which simplifies the encoding procedure and can 
be used for time-varying channel as well. The source 
uses superposition coding to 'route' part of its data 
through the relay. 

These aspects will be clear from the actual description of the 
scheme. Moreover we can show that 

Proposition 6. The LDA strategy achieves the gDoF upper 
bound in (|3]l. 

Proof: The proof can be found in Appendix [E] ■ 



A. A Motivating Example 

The LDA of the G-HD-RC in ([TJ is a deterministic channel 
with input-output relationship 



Yr 

Yd 



S"-^-X, (1-5,), 
S»-Ad^^ + S"-'^''^^, 5,, 



(29a) 
(29b) 



for some nonnegative integers /3sri Adj Ad^ where the inputs 
and outputs are vectors of length n := max{/3sr, /3sd, Ad} and 
S is the n X n shift matrix 110,1 . 

The capacity of a deterministic RC is given by the cut-set 
upper bound [TOl. For the LDA in (|29| the cut-set upper-bound 
evaluates to 
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^(HD) ^ f Ad + niax^g[o,i] mill |(1 - 9* (7)) log ^./.^^j + 6** (7) log ^(:^,7[/3sr - Ad] + } for /3sr > Ad, Ad > Ad 
[ Ad othei-wise. 

(30) 



max{i?}== max min X„ 5,; y^), /(X,; F,, F^l^r, 5r)| 

in{if(rd),ij(i;,rd|x,,5,)} 

min \^H{Yd\Sr),H{Yr, Yd\Xr, Sr)] + H{Sr) 



= max mm ■ 

PXs ,Xr.Sr 



< max 

< max min{7Ad + (1 - 7) max{Ad, Ad}, 7 max{Ad, Ar} + (1 - 7)Ad} + log(2) 

7G[0,1] 

= Ad + 7LDA[Ar - Ad]+ + l0g(2), (31) 



Theorem 2. The capacity of the deterministic HD RC in ( |29[ ) 

is given by at the top of this page where 9* (7) = 1 — 
max{l/L,7} and L := 2^^"^-^^-^^^ . 

Proof The proof can be found in Appendix IK] ■ 
Next we further upper bound the capacity in ( |30| because 
our goal is to get insights into asymptotically optimal strategies 
for the G-HD-RC. For the channel in (|29]) we have that ( |3T] ) 
at the top of this page holds, where 7£i3a. is the optimal 7 :— 
¥[Sr = 0] e [0,1] obtained by equating the two arguments 
within the min and is given by 

(A-d-fed 

7lda 



Remark 2. The HD optimal strategy in Figs. 2(a) and 2(b) 



should be compared with the FD optimal strategy in Fig. \2(c) 



(A, 




-Ad) + (Ar-Ad) 



for Ad > Ad, 
otherwise. 



Ar > A 



sd 



Next we show that the upper bound in ( |3T| l is achievable 
to within log(2) — 1 bit. This 1 bit represents the maximum 
amount of information I{Sr\Yj) that could be conveyed to 
the destination by a random switch at the relay. If we neglect 
the term log (2) we can achieve the upper bound in pT] ) 
with the scheme shown in Figs. 2(a) and 2(b) for the case 
min{Ar,Ad} > Ad, which is the case where the upper 
bound differs from direct transmission, i.e., X„ = 0. In 



Phase I/Fig. 2(a) the relay listens and the source sends hi (of 
length Ad bits) directly to the destination and 62 (of length 
Ar — Ad bits) to the relay; note that 62 is below the noise 
floor at the destination; the duration of Phase I is 7, hence 
the relay has accumulated 7(Ar — Ad) bits to forward to the 
destination. In Phase II/Fig. |2(b) the relay forwards the bits 
learnt in Phase I to the destination by 'repackaging' them into 
a (of length Ad — Ad bits); the source keeps sending a new 
hi (of length Ad bits) directly to the destination; note that a 
does not interfere at the destination with 62; the duration of 
Phase II is such that all the bits accumulated in Phase I can 
be delivered to the destination, that is 

7(Ar-Ad) = (i-7)(Ad-Ad), 

which gives precisely the optimal 7£qa- The total number of 
bits decoded at the destination is 

1 • Ad +7ldA • (Ar - Ad), 

which gives precisely the optimal gDoF for the half-duplex 
channel in ([3]). 



In Fig. 2(c) in a given time slot t, the source sends hi[t\ (of 
length Ad bits) directly to the destination and 62 + 1] (of 
length at most Ar — Ad bits) to the relay; the relay decodes 
both hi[t] and 62 [i + 1] and forwards 62 [i+ 1] in the next slot; 
in slot t the relay sends 62 [t\ ( of length at most Ad — Ad 
bits) to the destination; the number of bits the relay forwards 
must be the minimum among the number of bits the relay 
can decode (given by Ar — Ad) <^nd the number of bits that 
can be decoded at the destination without harming the direct 
transmission from the source (given by Ad — Ad)- Therefore, 
the total number of bits decoded at the destination is 

Ad + min{Ad - Ad, Ar - Ad}, 

which gives precisely the optimal gDoF for the full-duplex 
channel in Q 

Remark 3. Fig. ^compares the capacities of the FD and HD 
LDA channels; it also shows some achievable rates for the HD 
LDA channel. In particular, the capacity of the FD channel is 
given by ^ ( dotted black curve labeled "FD "), the capacity 
of the HD channel is given by ( |30[ ) ( solid black curve labeled 
"HD " obtained with the optimal pg in Appendix and its 
upper bound by (red curve labeled "HDlda upper"). For 
comparison we also show the performance when the source 
uses i.i.d. Bemoulli{l/2) bits and the relay uses one of the 
following strategies: i.i.d. Bernoulli(q) bits and random switch 
(blue curve labeled "HDiid q+rand" obtained by numerically 
optimizing q € [0,1]), i.i.d. Bernoulli{l/2) bits and random 
switch (green curve labeled "HDiid 1/2+rand" obtained with 
Pq ~ 1/L in Appendix and i.i.d. Bernoulli(l/2) bits and 
deterministic switch (magenta curve labeled "HDiid 1/2+det" 
and given by Ad + min{7[Ar - Ad]"*", (1 - 7)[Ad - Ad] + })- 
We can draw conclusions from Fig. [ij 

• With deterministic switch: i.i.d. Bernoulli{\/2) bits for 
the relay are optimal but this choice is quite far from 
capacity (magenta curve vs. solid black curve); this 
choice however is at most one bit from optimal (magenta 
curve vs. red curve). 

• With random switch: the optimal input distribution for the 
relay is not i.i.d. bits; i.i.d. inputs incurs a rate loss (blue 
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curve vs. solid black curve); if in addition we insist on 
i.i.d. Bernoulli{\/2) bits for the relay we incur a further 
loss (green curve vs. blue curve). 
This shows that for optimal performance the relay inputs are 
correlated and that random switch should be used. 

B. An achievable strategy inspired by the LDA 



We can mimic the LDA strategy in Section V-A for the G- 
HD-RC as follows. We assume S < C, otherwise we use direct 
transmission to achieve R — log(l + S). The transmission is 
divided into two phases: 

• Phase I of duration 7: the transmit signals are 

Xs [1] = VT^Xb, [1] + VsXb, , <5 ^ 

Xr[l] = 0. 



1 + 5' 



The relay applies successive decoding of [^j followed 
by Xb^ from 

Yr[l] = VC Vl-SXb.ll] + VC VSXb, + Zr[l], 

which is possible if (rates are normalized by the total 
duration of the two phases) 

<7log(l + C)-7log(^l + C-^ 
1 



S 



Rb, <7log l + C 



1 



(32) 



The destination decodes X^^yl] treating Xi,2 as noise from 

Yd[i] = Vs VT^Xh.ii] + Vs VsXb, + Zd[i], 

which is possible if 

(33) 



< 7 log (1 + 5) - 7 log (^1 + 5— 

Finally, since we assume S < C, Phase I is successful 
if ( (32] i and ([33| are satisfied. 

Phase II of duration 1 — 7: the transmit signals are 

Xs[2] = Xi,^[2] 

Xr[2] — Xb2 

The destination applies successive decoding of X^^ (by 
exploiting also the information about 62 that it gathered 
in the first phase) followed by -^^^[2] from 

Yd[2] = VSXi„i2] + c+^'VlXi,, + Zd[2], 
which is possible if 



(34) 
(35) 



Rb,< (1-7) log (i+i^)+7log(i+Y^ 

Rb,l2] < (l-7)log(l + 5). 

By imposing that the rate Ri,^ is the same in both phases, 
that is, that (|32]| and ([34]) are equal, we get that 7 should 
due chosen equal to 7* 



log 1 



I 

l+S 



log (1 + t:^) +l°g (1 + Ws) -log (1 + Ws) 



Note that 7* — > 7£]3a as SNR increases. Moreover 
we give here an explicit closed form expression for the 
optimal duration of the time the relay listens to the 
channel. 

The rate sent directly from the source to the destination, 
that is, the sum of (|33]) and dSSll, is 



R 



■bi 



■1] + Rb,[2] = l0g(l + 5) - 7* log (1 + 



e[oaog(2)] 



Therefore the total rate decoded at the destination through 
the two phases is 

Rb.ll] + Rb,m + = r(LDA-HD) ^ 

We notice that the rate expression for 7-(lda-hd) -^^ ^j]) , 
which was derived under the assumption C > 5, is valid 
for all C since for C < 5 it reduces to direct transmission 
from the source to the destination. 



VI. Single Relay Networks: Analytical Gaps 

In the previous sections we described upper and lower 
bounds to determine the gDoF of the G-HD-RC. Here we 
show that the same upper and lower bounds are to within a 
constant gap of one another thereby concluding the proof of 
Theorem [T] We consider both the case of random switch and 
of deterministic switch for the relay. 

Proposition 7. [PDF and random switch] PDF with random 
switch is optimal to within 1 bit. 

Proof The proof can be found in Appendix |G] ■ 

Proposition 8. [PDF and deterministic switch] PDF with 
deterministic switch is optimal to within 1 bit. 

Proof: The proof can be found in Appendix |H] ■ 

Proposition 9. [LDA (deterministic switch)] LDA is optimal 
to within 3 bits. 

Proof: The proof can be found in Appendix |l] ■ 

We conclude this section with a discussion on the gap that 
can be obtained with NNC. The NNC strategy is a network 
generalization of the CP. It has been proposed for general 
memoryless networks and it is optimal to within a constant gap 
for full-duplex multicast networks with an arbitrary number 
of relays, where the gap grows linearly with the number of 
relays ifTTI . In the case with only one relay, NNC reduces 
to the classical CP Remark 18.6] and represents a good 
alternative to the PDP especially in the case when the link 
between the source and the relay is weaker than the direct 
link. The NNC rate is presented in Appendix |F] By using 
Remark |5] in Appendix |F| we have 

Proposition 10. [NNC and deterministic switch] NNC with 
deterministic switch is optimal to within L61 bits. 

Proof: The proof can be found in Appendix |j] ■ 



g 



VII. Single Relay Networks: Numerical Gaps 

In this section we show that the gap results obtained in 
Section [Vl] are pessimistic and are due to crude bounding in 
both the upper and lower bounds, which was necessary in order 
to obtain rate expressions that could be handled analytically. In 
order to illustrate our point, we first consider a relay network 
without the source-destination link, that is, with 5 = 0, in 
Section IVII-AI and then we show that the same observations 



The LDA achievable rate in ( fTT) with S" = is given by 



^(LDA-HD)|^^^ 



are valid for any network in Section VII-B 



A. Single Relay Networks without a Source-Destination Link, 
a.k.a. Diamond Networks with One Relay 

c) Upper Bound: We start by showing that the (upper 
bound on the) cut-set upper bound in (|6]) can be improved 
upon. Note that we were not able to evaluate the actual cut- 
set upper bound in (|5]l so we further bounded it as in (|6|, 
which for 5 = reduces to 



max min j'H(7) + (l — 7) log ( IH — — 
7e[oa] [ V 



7log 1 



C 

7 

The capacity of the G-FD-RC for 5 = is known exactly and 
is given by the cut-set upper bound 

C(™)|s=o=log(l + min{C,/}). 

is a trivial upper bound for the capacity of the G-HD- 
RC. Now we show that our upper bound r^'-'^^^'-*) [5=0 can 
be larger than C'^™^s=a- For the case C = 15/2 > I = 3/2 
we have 

r(^^-™) |s=o > minj^ Q + ^ log (1 + 2/) , ^log (1 + 26) | 

= log(4)>C7(™)|s=o-log(2.5). 

The reason why the capacity of the FD channel can be 
smaller than our upper bound r^^^^^^^js^o is the crude 
bound I{Sr;Yd) < H{Sr) = ^(7). As mentioned earlier, we 
needed this bound in order to have an analytical expression for 
the upper bound. Actually for 5 = the cut-set upper bound 
in (|5]l is tight, as we show next. 

d) Exact capacity with PDF: 

Theorem 3. In absence of direct link between the source and 
the destination PDF with random switch achieves the cut-set 
upper bound. 

Proof: With 5 = 0, the cut-set upper bound in (|5]l and the 
PDF lower bound in ([8| are the same (see also Appendix |G] 
with 5 = 0). ■ 

e) Improved gap for the LDA Lower Bound: Despite 
knowing the capacity expression for 5 = 0, its actual eval- 
uation is elusive as it is not clear what the optimal input 
distribution ^ g in (|5]l is. For this reason we next 
specialized the LDA strategy to the case 5 = and evaluate 
its gap from the (upper bound on the) cut-set upper bound 
in (|6]l. 



max min{7log (1 + C) 

76[oa] 

(l-7)log(l + /)}. 

and its gap from the outer bound can be reduced from 3 bits 
to about 1.5 bits since 

GAP < r(cs-HD)|^^^ _ ^(lda-hd)|^^^ 

< max |7log (1 + — | - 7 log (1 + C) , 
7e[oa] [ V 7/ 

H(7) + (1 - 7) log (1 + - (1 -t) log (1+/) 

< max j7log(- ) ,-H(7) + (l-7)log(— 
7e[oa] I V7/ Vl -7 

= max |-H(7) + (l-7)log f ^— ) 1 = 1.5112 bits. 
7e[04] t Vl~7/J 

Note that the actual gap is even less than 1.5 bits. 
By numerically evaluating the difference between 

min{C(FD),r(cs-HD)||^^^ ^(lda-hd) j^^^ ^und 

that the gap is at most 1.11 bits. 

f) Numerical gaps with deterministic switch: Similarly, 
by numerical evaluations one can find that the PDF strat- 
egy with deterministic switch in Remark |4]- Appendix |C] and 
the NNC strategy with deterministic switch in Remark |5]- 
Appendix [F] are to within 0.80 bits and 1.01 bits, respectively, 
of the improved upper bound. Notice that in these cases there 
is no information conveyed by the relay to the destination 
through the switch. Further reductions in the gap with random 
switch are discussed next for a general network. 

Fig. |4] shows different upper an lower bounds for the G- 
HD-RC for 5 = 0, C = 15, / = 3 vs 7 = V[Sr = 0]. 
We see that the cut-set upper bound exceeds the capacity 
of the G-FD-RC (maximum of the solid black curve vs. 
dashed black curve). Different achievable strategies are also 
shown, whose order from the most performing to the least 
performing is: PDF with random switch (red curve, 1.913 
bits/ch.use), PDF with deterministic switch (blue curve, 1 .702 
bits/ch.use), NNC with random switch (cyan curve, 1.446 
bits/ch.use), NNC with deterministic switch (magenta curve, 
1.402 bits/ch.use), and LDA (green curve, 1.333 bits/ch.use). 
In this particular setting, the maximum rate using the NNC 
strategy with random switch (cyan curve, 1.446 bits/ch.use) 
is achieved for V[Q = 0, 5^ = 0] = 0,P[Q = 0, 5,. = 1] = 
0.33, P[Q = 1,5^ = 0] = 0.45, P[g = 1,5^ = 1] = 0.22. 
This is due to the absence of the direct link (5 = 0) between 
the source and the destination. Actually, since the source can 
communicate with the destination only through the relay, it 
is necessary a coordination between the transmissions of the 
source and those of the relay. This coordination is possible 
thanks to the time-sharing random variable Q, i.e. when Q = 
the source stays silent, while when Q = 1 the source transmits. 

B. Single Relay Network with a Source-Destination Link 



Although the considerations in Section VII-A were for a 
relay channel without a source-destination link, they are vaUd 
in general. 
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Fig. |5] and Fig. |6] show the rates achieved by using the 
different achievable schemes presented in the previous sections 
with S' > 0. In Fig. |5] the channel conditions are such that 
the PDF strategy outperforms the NNC, while in Fig. |6] the 
opposite holds. In Fig. |5] the PDF strategy with random switch 
(red curve, 11.66 bits/ch.use) outperforms both the NNC 
with random switch (cyan curve, 11.11 bits/ch.use) and the 
PDF with deterministic switch (blue curve, 11.4 bits/ch.use); 
then the PDF with deterministic switch outperforms the NNC 
with deterministic switch (magenta curve, 10.94 bits/ch.use), 
which is also encompassed by the NNC with random switch. 
Differently from the case without direct link, we observe that 
the maximum NNC rates both in Fig. [5] and in Fig. |6] are 
achieved with the choice Q = %, i.e. the time-sharing random 
variable Q is a constant. This is due to the fact that the 
source is always heard by the destination even when the relay 
transmits so there is no need for the source to remain silent 
when the relay sends. 

In Fig. [7] we consider the case of deterministic switch. Fig. [7] 
shows, as a function of SNR for /3sd = 1, (Ad, Ar) G [0, 2.4], 
the maximum gap between the cut-set upper bound 7'(cs-hd) 
in (|6| and the following lower bounds with deterministic 
switch: the PDF lower bound obtained from 7-(pdf-hd) -^^ ^ 

fPDF) I 1^ — 

with /g — 0, the NNC lower bound in Remark pi in 
Appendix |F[ and the LDA lower bound in ( fTTj i. From Fig. |7] 
we observe that the maximum gap with PDF is of 1 bit as 
in Proposition [8] but with NNC is around 1.52 bits and with 
LDA is around 1.59 bits, which are lower than the analytical 
gaps we found in Propositions 10 and [9] respectively. 



The lower bounds can be improved upon by considering 
that information can be transmitted through a random switch 
for the relay. However, this improvement depends on the 
channel gains. If the information cannot be routed through 
the relay because min{C, /} < S, then the system cannot 
exploit the randomness of the switch, and so /^^^ = and 
jNNC — are approximately optimal (in this case the relay 
can remain silent). For this reason the maximum numerical 
gap obtained with a random switch coincides with the one 
obtained with a deterministic switch, as there are channel 
conditions for which random switch is not necessary. This 
behavior for the PDF strategy is represented in Fig. [8] In 
this figure we numerically evaluate the difference between the 
analytical gap, i.e., the one computed with 1^^^ = 0, and 
the numerical one, i.e., computed with /q"^* (actual value of 
/^°F), at a fix SNR = 20dB and by varying (/3rd,/3sr). We 
observe that when the information cannot be conveyed through 
the relay, i.e., min {Ad, Ar} < 1, then Iq^^ = is optimal, 
since the information only flows through the direct link. In 
Fig. |9] the channel channel gains are set such that the use of 
the relay increases the gDoF of the channel (Ad — 1 and 
(Ad,Ar) G [1,2.4]). Here the relay uses PDF. We observe 
that we have a further improvement in terms of gap by using 
a random switch (blue curve) instead of using a deterministic 
switch (red curve). We notice that at high SNR, where the 
gap is maximum, this improvement is around 0.1 bits. As 
mentioned earlier, the rate advantage of random switch over 
deterministic switch depends on the channel gains. 



VIII. Networks with Multiple Relays 

In this section we extend our gDoF and gap results to 
general HD-MRC. Similarly to the the full-duplex case [11], 
our main result is that NNC is optimal to within a constant gap 
for the HD-MRC where the gap is a function of the number 
of relays. 



A. Network Model 

In this model we have K > nodes, i.e., one transmit- 
ter/node 1, one receiver/node K and K — 2 relays with indices 
2, . . . ,_fir — 1. Each node is equipped with a single antenna 
and is subject to an average power constraint, which we set 
to 1 without loss of generality, i.e. E[|Xfcp] < 1 G M+, 
fee [\ : K — \\. The system is described by the input/output 
relationship 

Y = (I - S)HS X + Z (36a) 
Z = [Zi,...,Zk]^-A/'(0,I) (36b) 

Y = [yi,...,y;f]^ec^ (36c) 

X= [Xi,...,Xk]^GC^ : 

E[|^fcP] < 1 for fc e [1 : i-T - 1], (36d) 
H e C-^''^ (36e) 
S = diag{S'i, . . . , Sk} ■ 

Si ==1, Ske {0, 1} for fc e [2 : - 1], Sk = 0, (36f) 



where the vector S represents the state of the nodes, either 
receive (S — 0) or transmit (S ~ 1). The channel matrix 
H e C^^^ is constant and therefore known to all terminals. 
The entry hij with G [1 : K]'^ represents the channel 

between source j and destination i. Without loss of generality 
we assume that the noises Z are zero-mean jointly Gaussian 
and with unit power. Furthermore we assume that the noises 
are iid JV{0, 1), this is however not without loss of general- 
ity |19|. 

The capacity of the channel described in ( |36| ) is not known 
in general. Here we show that a scheme based on the NNC 
strategy achieves the capacity within a constant gap for any 
number of relays and for any choice of channel parameters. 
As for the FD case, the gap is found to be a function of the 
number of relays. We propose a different and simpler bounding 
technique than that of 1.1 1 J , which might overestimate the 
actual gap between inner and outer bound. 

B. Capacity to within a constant gap 
Our main result is 

Theorem 4. The cut-set upper bound for the half-duplex multi- 
relay network is achievable to within 

GAP < max { minil + £, /sT - 1 - ^} log (1 + ^) 

£e [0:^-2] 

+ mm{l + UJ. + K -I}} (37) 

bits per channel use. 

Proof: The proof can be found in Appendix |L] ■ 
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For high value of K, i.e. K >> 1, the optimal value of £ channel gains as 



in 



IS I = 



A'-2 



and the corresponding limit gap becomes 



GAP 



K 



\og(4K). 



(38) 



Fig. 10 shows the trend of the gap in ( |37| i as a function of the 
number of nodes K (blue curve). The red curve represents the 
limit behaviour of the gap in ( |37] i. 

Note that by applying the above gap-result to the case K = 
3 we obtain GAP < 4 bits which is a much larger gap than 
1.61 bits we found for the case of one relay. 

A smaller gap may be obtained by computing tighter 
bounds. This can be accomplished by several means. For 
example: 

• By using more complex and sophisticated bounding tech- 
niques: in 1 11 1 an upper bound on the water-filling power 
allocation for a general MIMO channel was derived. This 
more involved upper bound could be used here to obtain a 
smaller gap. We note that our bounding technique applied 
to the FD-MRC gives a gap of order f log(2i^), which 
is larger that 0.63/^ found in fTTl. 

• By using an achievable strategy based on PDF, which in a 
single relay case gives a smaller gap than NNC. However, 
PDF seems not to be easily extended to networks with an 
arbitrary number of relays, which is the main motivation 
we consider NNC here. 

• By deriving tighter bounds on specific network topolo- 
gies: in ll20j it is found that for a FD diamond network 
with K relays the gap is of the order \og{K), rather than 
linear in K fTTl. Moreover, for a symmetric FD diamond 
network with K relays the gap does not depend on the 
number of relays and is upper bounded by 2 bits. The 
key difference between a general multi-relay network and 
a diamond network is that for each subset A we have 
Rank[H^ = 2, i.e., the rank of a generic channel sub- 
matrix does no longer depend on the cardinality of A. 
Based on this observation we have 

Proposition 11. The cut-set upper bound for the gaussian 
half-duplex diamond network with {K — 2) relays is 
achievable to within 

GAP <{K -2) log(2) + 41og(if) + 21og(e/2). 

Moreover, if the conjecture in H17\l holds then the gap 
above could be decreased to 

GAP < 51og(/0 + 21og(e/2). 

Proof: The proof can be found in Appendix [M] ■ 



C. Example: Fully connected network with K = 4 

To gain insights into how relays are best utilized, we 
consider a network with two relays, i.e. K ^ 4 nodes. In 
particular we highlight under which channel conditions the 
gDoF performance is enhanced by exploiting both relays rather 
than using only the best one. We also compare the loss that 
incurs by using HD with respect to FD. Let parameterize the 



log(|/».,f) 



log(SNR) 



(^J)6[l:4]^ 



/ * 

asi 

V 1 



132 

aid 



Pi 
a2d 



) 



where * denotes entries that do not matter; this is so because 
the source node never listens to the channel (first row), the 
destination node never transmits (last column), and a relay 
can remove the 'self-interference' (main diagonal). Note that 
the direct link from the source to the destination has gain SNR 
and all other channel gains are expressed with reference to it. 

g) Full Duplex: In the following the set A indicates the 
relays that lie on the source/node 1 side of the cut. The cut-set 
upper bound for the FD-MRC LlJj gives 



A = 



A ^{2} 



A ^{3} 



^ = {2,3} 



h ■.= I{X,;Yi,Y^,Y3\X2,X3) 

<l0g(l + |/l41P + |/l21P + |/i3in 

/2 ■.= I{Xi,X2;Ya,Y3\X3) 

<log|l2+HiHf|+log(2) 

= log(l + A)+21og(2) 
h ■.^I{Xi,X3;Yi,Y2\X2) 

<log|l2+H2Hf|+log(2) 

= log(l + i?) + 21og(2) 
h :-/(Xi,X2,X3;r4|0) 

<l0g(l + (|/l4l| + |/l42| 



where 

A- 
B 



i-2 = 



hAiV- 

/l41P 

1 0' 
1 



'•421 
|/i43|' 



Hi 



\h2i\^ + \ 

h^l /l32 
^41 



\h,3\y 



/132I + 1^31^42 - ^32 ^41 1 
/l23p + 1^21^43 - ft-23/l41P 
/^21 h23 



and H2 = 



/141 /14 



Thus for the full-duplex case, the cut-set bound is given by 

which at high SNR gives the following gDoF (achievable to 
within a constant gap 111]) 



h(fd) 



lim 



^(CS-FD) 



SNR^+00 log(l + SNR) 
= min I max{l,Q;si, Q!s2} , max {0^2 + aid, /32 + 1} , 
max{asi + a2d, Pi + 1} ,max{l, aid, a2d} }• 

In the following we are interested in identifying cases for 
which 



d]f=4 > dK=4,best relay := Hiax 1 1 , min{a,i , aid} , 

min{as2,Q!2d}}, 



where d 



(FD) 
Ar=4,be.st relay 



is the gDoF that can be obtained by 



selecting the relay which achieves the highest gDoF in (HI) 
while leaving the other silent. From the expression for 
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we immediately see that in order to have d^^| > 1, i.e., better 
than direct transmission with the two relays silent, we need 

max{Q;si, q;s2} > 1, 
max {as2 + aid - ^,^2} > 0, 
max{asi + a2d - 1, /3i} > 0, 

max{aid,a2<i} > 1- 

We therefore distinguish the following cases: 

1) if asi > max{l,as2} and aid > max{l,a2rf} then 



OK=i = min{Q;si,aid} 



~ '-'_fS'=4,best relay ^ ^' 



therefore this case is not interesting; 

2) if asi > max{l,as2} and a2d > max{l, 

^K=l = niin|asi,max{Q!s2 + aid, (^2 
Next consider the following sub-cases: 



aid} then 
+ l},a2<i|- 



< 1, then d 



(FD) 



< d 



(FD) 



a) If aid < 1 or a,2 - ^, UK=4,best relay ^K=4' 

hence in this case using both relays can give an 
unbounded rate improvement over using the best 
relay if max {as2 + aid ~ 1, /32} > 0; 

b) If aid > 1 and 0^2 

max{ min{asi,aid},min{as2,a2d}} > 1; we 
distinguish 

i) If a^i < aid or a2d < then 



> 1, then d^f] best relay 



j(FD) 
if =4, best relay 



l{asl,a2d} = d 



(FD) 



therefore this case is not interesting; 
ii) If asi > aid and a2d > C(s2 then 

j(FD) r 1 ^ j(FD) 

dk-=4,best relay = max{aid,as2} < d)^^^, 
therefore we have a strict improvement by 
using both relays over using only the best relay. 
3) Considering cases similar to the two above but with the 

role of the relays swapped complete the list of possible 

cases. 



An example of network satisfying the conditions in item 2a 
or item 2(b)ii is given in Fig. 11 where the numerical value 



on a link represents the SNR exponent on the corresponding 
link. 

h) Half Duplex: With HD, each term in the min-function 
defining the cut-set bound is a convex combination of 2^~^ 
terms, one for each possible state of the network. Formally, 
these terms are as for the FD case but where one replaces the 
SNR exponent of the link from node j to node i, given by 
aij, with (1 — Si)aijSj, where {Si,Sj) £ {0,1}^ indicates 
the state of the nodes. Because no additional insight is gained 
from these terms, we do not report them here. At high SNR 
the outer bound gives the following gDoF 



j(HD) 

d)^^4 = max mm 



[XooDf^+XoiD^i^+XioD^i^+XiiD^^^ , 



^00^2 +^oi-t^2 +MqJ-^ 



XooD''" +X01D''' +XioD'^' +X11D, 



.(0) 
^3 

^4 



XooDl^'+XoiDi'hXioD'/'+XiiD, 



2 +>^iiD^2 

(2 
3 

(2) 

4 



(3) 



(3) 
3 

(3) 
4 



where the maximization is over Xij, € {0, 1}^ represent- 
ing the fraction of time node 2 is in state 82^1 and node 3 



is in state 5*3 = j 



D 



(0 



(1 



D 



(2 



D. 



(3 



D, 
D 
D 
D 



such that Aoo + Aoi + Aio + An = 1 and 
:= max{l,asi,as2}, 



of'' := max{l,asi} 



D^^^ := max{l,as2} : 



= 1, 

max {as2 + aid, ^2 + ^ , 



D 



(2) 



max {1, aid} : 



max {ttsi + a2d, /3i + 1} , 



D 



(1) 



max{l,a2d} : 



max{l,aid,a2d} ■ 



An analytical closed form solution for the optimal {A^} is 
complex to find for general channel gains. However, numer- 
ically it is a question of solving a linear program, for which 
efficient numerical routines exist. We remark that this linear 
program can be thought as the high SNR solution of the 
iterative algorithm proposed in |18|. 

Table |l] shows the gDoF corresponding to the four cases 
listed in the FD case where using both relays strictly im- 
proves over exploiting only the best relay. We denote with 
diJ=4 best relay gDoF that is obtained when the two relays 
work in HD by selecting the relay which achieves the highest 
gDoF while leaving the other silent, whose closed form 
solution is given in ([3]l. From Table I!] we notice that in each 
case also for the HD case we have d^^} > d^^^l best relay ^ 
for the FD case. Furthermore, as expected, d^^| > d^^^. 

IX. Conclusions 

In this work we considered the HD-RC, a network where 
the source communicates with the destination through a relay 
node. The relay node works in half-duplex, in the sense 
that it transmits and receives in different time slots. This 
scenario, with the half-duplex assumption, represents a more 
practically relevant model compared to the full-duplex case. 
In particular, it is applicable in practical relaying architectures 
for 4G cellular networks. 

We derived in a close form expression the generalized 
Degrees-of-Freedom of this system and we show that three 
schemes achieve the capacity to within a constant gap. The 
first scheme turns out to be simple both in the encoding and 
decoding phases and it is inspired by the Linear Deterministic 
Approximation of the Gaussian noise at high SNR; the second 
and the third schemes are based on the Partial-Decode-and- 
Forward and Noisy-Network-Coding strategies. 

All these schemes consider both deterministic and random 
switch at the relay. In the first case the switch is known by 
all the terminals in the network, while in the second case 
the randomness that lies into the switch can be exploited to 
send more information and so to achieve higher rates and 
consequently lower gaps. We show that random switch is 
optimal and in some cases, as in the diamond network with 
only one relay, achieves the exact capacity. 

Finally we extend our results to a multi-relay channel and 
we show that a strategy based on the NNC with random switch 
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TABLE I: gDoF when using both relays is better than using only the best one 



Channel Parameters 


Full-duplex 
best relay 


Full-duplex 
both relays 


Half-duplex 
best relay 


Half-duplex 
both relays 


(2.5,1.4,0.5,1.8,0.6,0.8) 


1.4 


1.8 


1.267 


1.4235 


(2.5,0.3,0.7,1.3,0.4,0.8) 


1.0 


1.3 


1.000 


1.2182 


(1.8,1.2,1.3,2.0,0.7,1.2) 


1.3 


1.8 


1.218 


1.5808 


(1.7,1.1,1.2,1.4,0.4,1.5) 


1.2 


1.4 


1.156 


1.3604 



at each relay, achieves the capacity to within a constant gap. 
We also prove that this gap may be even decreased in more 
structured settings as, for example, the diamond network where 
there are no source-destination and relay-relay links. 
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Appendix A 
Proof of Proposition[T] 

Proof: An outer bound to the capacity of the memoryless 
RC is given by the cut-set outer bound |4, Thm.16.1] that 
specialized to our G-HD-RC channel gives 



-RC) 



^(HD 



< max min 



in {l{Xs, [Xr ,Sr];Yd), I{X, ; F, , Yd \ [Xr , 5,-] ) } 



max mm 

PXs,Xr,S. 



< max min 



< max min 

^. ^(CS-HD 



m{l{Sr;Yd)+I{Xs,Xr;Yd\Sr), 

li^s'iYr, Yd\Xr, Sr)^ 
{H{Sr)+I{X,,Xr;Yd\Sr), 
I{Xs;Yr,Yd\Xr,Sr)} 

{n{-f) + + (1 - -f)l2, ih + ii 



(39a) 
(39b) 

(39c) 

7^4 } 

(39d) 



where the different steps follow since: 

• We indicate the (unknown) distribution that maxi- 
mizes P9a[ ) as X s order to get the bound in (|5]l. 

• In order to obtain the bound in ([39c| we used the fact 
that for a discrete binary-valued random variable Sr we 
have 

I{Sr;Yd) = H{Sr) - H{Sr\Yd) < H{Sr) = Hin) 

for some 7 :— P[S',. = 0] G [0, 1] that represents the 
fraction of time the relay listens and where 7^(7) is the 
binary entropy function in (T2|. In (39d| the maximization 
is over the set defined by ( T3]l-([T5]l and is obtained as an 
application of the 'Gaussian maximizes entropy' principle 
as follows. Given any input distribution Pxs,x^.S-r^ '^he 
covariance matrix of [Xg, Xj.) conditioned on Sr can be 
written as 



Gov 



Ps,i Oil ^Ps,lPr,l 

\/ Ps,ePr,£ Pr,£ 



with |a^| < 1 for some {Ps,o, Ps,i, Pr.Oi Pr.i) G sat- 
isfying the average power constraint in ( [T5] l. Then, a zero- 
mean jointly Gaussian input with the above covariance 
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matrix maximizes the different mutual information terms 
in ( |39c| ). In particular 

IiXs,Xr; Yd\Sr = 0)< log (1 + SP,,o) h, 

I{Xs,Xr;Ya\Sr = 1) 

<l0g {l + SPs,l+IPr,l+2\aiWSPs,l IPr,l) =-h, 
I{X,-Yr,Yd\Xr,Sr=Q) 

<log(l + (C + ^)(l-|aonPs,o) 
<log(l + (C + 5)P,,o) =:/3, 

I{X,;Yr, Yd\Xr, 5r = 1) < log (l + S{1 - |ai|')Ps,l) 

=:/4, 

as defined in ([T6]l-(|l9]l thereby proving the upper bound 
in (|6]l, which is the same as r^^^^^^' in P9d[ ). This 
shows the bound in (|6]l. 

In order to get to (|7]i from (|6]) we let the channel gains be 
parameterized as in (|2]i. The average power constraints at 
the source and at the relay given in ([Tsj can be expressed 
as follows. Since the source transmits in both phases we 
define for some /3 € [0,1] 



Psfi — 



Ps,i - 



7 

1-/3 
1-7' 



On the other hand, the relay transmission only affects 
the destination output for a fraction (1 — 7) of the time, 
i.e., when 5*^ = 1, hence the relay must exploit all its 
available power when S'r = 1; we therefore let 

Prfi = 0, 

1 

Pr,l = 



1-7 



With this, the cut-set upper bound r.(CS-HD) p9d[ ) 
can be rewritten as ( |40| l at the top of next page where we 
defined bi and 62 as in (|20|-(|2T]i, namely 

\og(i + {Vi + Vs)A 

bi ■■= 7T— > 1 since / > 0, 

log (1 + S) 

\og{l + C + S) 



log (1 + 5) 



> 1 since C > 0. 



Note that the optimal 7 is found by equating the two 
arguments of the maxmin and is given by 



7cs := 



(bi 1) 



(fol-l) + (fo2-l)' 

This proves the upper bound in (|7|. 



Appendix B 
Proof of Proposition[2] 

Proof: The upper bound in (|7| implies 

^(HD-RC) 

log (1 + 5) / , (6i-l)(62-l) 



< lim 



1 



SNR^+00 log (1 + SNR) V (bi - 1) + (^2 - 1) 

[/3rd//3sd - 1]+ [Psr/Psd - 1] + 
[/3rd//3sd - 1]+ + [Psr/Psd - 1] + 
[/3rd - /3sd]+ Wsr - /3sd] + 



= /3sd [I + 
= /3sd + 



[/3rd - ^sd]+ + Wsr - /3sd] 



+ ' 



since bi max{/3sd, /3rd}//3sd and 62 max{/3sd, Ar}/;Ssd 
at high SNR, which is equivalent to the right hand side of ^ 
after straightforward manipulations. ■ 



Appendix C 
Proof of Proposition[3] 

The largest achievable rate for the memoryless relay channel 
is the combination of PDF and CF/NNC proposed in the 
seminal work of Cover and ElGamal |3|. Here we use PDF. 
The bound in (|9]l follows since: 

Proof: The PDF scheme in f4^, Thm. 16.3] adapted to the 
HD model gives the following rate lower bound 

^(HD-RC) 

> max mm\liSr]Ya)+IiXs,Xr;Yd\Sr): 

I{U- Yr\Xr, Sr) + /(X.; Yd\U , Xr, Sr)] 

> maxmin {4''°^' + jh + (1 - 7)^,7^ + (1 " 7)^3} 
= in (|9J, 

where for the last inequality we let 7 := ¥[Sr = 0] e [0, 1] 
be the fraction of time the relay listens and, conditioned on 
Sr = £ E {0, 1}, we consider the following jointly Gaussian 
input 



0, 



■ 1 

.Pr\l 



Pt\l Pr\e 
1 at 

1 



a 



1 Pt\e Pr\i 
Pt\e 1 



1 



>- 0. 



In particular, we use specific values for the parameters 

{Pt\i,Pr\i,ai}i<z{o^i}, namely 

Zai+0 = O, (41a) 
ao = and either |pt|o|^ = 1 - |pr|o|^ = 

or \prio\^ = 1 - |pt|o|2 = 0, (41b) 

Pt\i = al, Pr\i = 1- (41c) 
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^(CS-HD) ^ 



(7,|ai|,/3)e[o,i]-^ 1^ V 7/ \l-7 1-7 yi-71-7 



+ 7log ( 1 + ^ + ^"l + (1 - 7) log f 1 + (1 - 



7 7/ V 1 — 7 

2N 



< max niin 

7e[o,i] 



{H(7)+7log(l + ^)+(l-7)log 
+ 7 log 1 + — + - + (1 - 7) log 1 + 



11) V 1 — 7 



max mm 

75 [0,1] 



in \m{-i) + 7 log (7 + 5*) + (1 - 7) log (^1 - 7 + (\// + TS^) 
H(7) + 7 log (7 + C + 5) + (1 - 7) log (1 - 7 + S)} 



< 2 log(2) + max^ min {7 log (1 + S*) + (1 - 7) log (^1 + (\/7 + ^ j , 7 log (1 + C + 5) + (1 - 7) log (1 + S*)} 

= 21og(2)+log(l + 5) max min {7 + (1 - 7)61, 763 + (1 - 7)} 

7e[o,i] 



21og(2)+log(l + 5) 1+ max min{(l - 7)(fei - 1),7(62 - 1)} 
V 7e[o,i] 



With these definitions, the mutual information terms we let U = Xr, we have 



(PDF) 







1 ^ o ; ■ ■ • ! 



h in ^ 



I{Xs,Xr;Yd\Sr = 0) = log(l + SPs,o) h; 

I{^s,Xr;Yd\Sr = 1) 



= log [1 + SPs^i + IPr^ + 2|ai|^5P,,i IPr,i) =:Ie, 

(note I5 = Ii and /g = /2 because of the assumption in ( |41a| l); 
next, by using the assumption in ( |41b[ ), that is, in state 5,. = 
the inputs Xg and Xr are independent, and that either U — Xg 
or U = Xr, we have: if U = Xs independent of Xr 

IiU;Yr\Xr,Sr^O)+I{Xg;Yd\U,Xr,Sr^O) = 
= I{Xg; VCXg + Zr\Xr, Sr = 0) 
+ I{Xs;VSXs +Zd\Xs,Xr,Sr^O) 

= log(l + CP,,o)+0, 
and if U — Xr independent of Xg 

IiU;Yr\Xr,Sr = 0) + /(X, ; Frf | [/, X„ 5,, - 0) = 

= I{Xr; VCXg + Zr\Xr,Sr = 0) 
+ I{Xs;^/SXs + Zd\Xr, Sr = 0) 

= o + iog(i + 5P,,o); 

therefore under the assumption in ( |41b[ ) we have 

I{U;Yr\Xr,Sr^O)+I{Xg;Yd\U,Xr,Sr=0) 

= log (1 + max{C, 5}P,,o) 
next, by using the assumption in ( |41c[ ), that is, in state Sr = 1 



/([/; Yr\Xr, Sr = l)+ I{Xg;Yd\U, Xr, Sr = 1) 
= I{Xr; Zr\Xr, Sr = 1) + I{X,; VSX, + Zd\Xr, Sr = 1) 
= + I{Xs; ^Xs + Zd\Xr, Sr = 1) 

= log(l + ^(l-|ainP,,i) =:Is, 
(note Ij < I J, and /g = I 4); finally 



l{Sr;Yd) 



E 



loe 



1 

fYAYd) 



' [7 log(wo) + (1 - 7) log(wi) + log(7re)] 



(PDF) 



where /y^ ( ) is the density of the destination output 1^, which 
is a mixture of (proper complex) Gaussian random variables, 
i.e., 

fvAt) = — exp(-|i| V«o) + — cxp(-|t|V«i), t G C, 
■kvq nvi 

vo Var[rd|^, = 0] =cxp(/5), 
VI Var[rd|5, = 1] = exp(/6). 

Note that 4^^^^ = < H{Sr) = ^(7). This proves 

the lower bound in (|9|. 

Next we show how to further lower bound the rate in (|9]) 
to obtain the rate expression in ( [TOj i. With the same parame- 
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(PDF-HD) • f r(PDF) , 1 A , 1^^ 

r*- = max mm /A '+ 7 log IH 

7e[o,i],|a|<i,/3e[o,i] I ^ V 7 



/ 5(1-/3) / , /S-fl-/?) / 

+ (l-7)log 1 + ^ ^ + - +2|ai ' ^ '^^ 



1 — 7 1 — 7 V 1 — 7I — 7 

(] 

1-7 



7log( 1+ ^max{C/3,5/?}') +(l-7)logfl + (l-|«n^^^ 



> max mill < 

7G[o,i],;3e[oa] 



7 log ( 1 + ^ max {/3C, pS}] + (1 - 7) log (l + 



7 / V (1-7) 

> max min {7 log (1 + 5) + (1 - 7) log (1 + S* + /) , 7 log (1 + max {C, S*}) + (1 - 7) log (1 + S*)} 

7G[0,1] 

= log(l + S) max min {7 + (1 - 7)01,702 + (1 - 7)} 
7e[o,i] 



log (1 + 5) 1+ max min {(1 - 7)(ci - 1), 7(02 - 1)} 
V 7e[o,i] 

log(l + 5) (l+ ^''-l^^'r'l )^ (42) 



(ci-l) + (c2-l);' 



terization of the powers as in Appendix |Aj namely 



Ps,0 — 
^..1 = 

Pr.l = 



7 

1-/3 
1 
0, 

1 

1^ 



7 



we have ( [42| i at the top of this page where we defined ci and 
C2 as in (|27]i-(|28]l, namely 



Appendix D 
Proof of Proposition^ 

Proof: The lower bound in ( [TO] i implies 



d(HD 



-RC 



-s^+00 log(l+SNR) 



(C1-1)(C2-1) 
(Cl - 1) + (C2 - 1 

[/3rd//3sd - 1]+ [Ar/Ad - 1]^ 



[/3rd//3sd - 1]- 
[Ad - Ad]+ [A: 



- [Ar//3sd - 1]- 



sd 



[Ad - Ad]+ + [/3sr - A 



sd 



since ci ^- max{/3sd, Adj/Ad and C2 ^ max{/?sd, Arj/Ad 
at high SNR, which is equivalent to the right hand side of ([3]) 
after straightforward manipulations. ■ 



\og{l + I + S) . 

Cl := — > 1 since / > 0, 

log(l + 5) 

log(l + max{C,g}) ^ ^ . 

C2 := 7- > 1 since C > 0. 

log (1 + 5) 

Notice that Ci < bi,i = 1,2, where bi,i = 1,2, are defined 
in pOl)-pT|). The optimal 7, indicated by 7pDp is given by 



7pdf := 



(ci 1) 



(Ci - 1) + (C2 - 1) 



e [0,1]. 



Remark 4. A further lower bound on the PDF rate 
j,(PDF-HD) jjgi obtained by trivially lower bounding 

fPDF^ ' — ' 

/g > 0, which corresponds to a fixed transmit/receive 
schedule for the relay. 



Appendix E 
Proof of Proposition^ 

Proof: The rate in ( fTT) can be further lower bounded as 



^(LDA-HD) > _ iQg(2) + log (1 + 5) ( 1 



(C3-1)(C4-1) 
(C3-1) + (C4-1) 

, l0g(l + / + S) J , l0g(l + C+5) rp. 

where cg := ci = f,;,(i+g) and C4 := 62 = rog(i+g) ■ The 
rate above implies 



d> lim ^"g + ^) 
- s4?oolog(l + SNR) 



= /3sd I 1 
-/3sd 



(C3-1)(C4-1) 
(C3 - 1) + (C4 - 1) 
[/3rd//3sd - 1]+ [/3sr//3sd - 1]^ 



[Ad/Ad - 1]+ + [/3sr//3sd 
[Ad - /3sd]+ [/3sr - /3sd] + 



11 + 



[/3rd 



/3sd]+ + [/3sr - /3sd] 



since C3 ^- max{/3sd, /3rd}//3sd and C4 max{Ad, /3sr}//3sd 
at high SNR, which is equivalent to the right hand side of ([3]l 
after straightforward manipulations. ■ 
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^(HD-RC)> mm\liX,;Yr,Yd\[Xr,Sr],Q), 
I{X, , [Xr , Sr] ; Frf I Q) - /(n ; f; |X„ [X, ,Sr],Yd,Q)} 



= max min < I{Xs;Yr,Yii\Q, Sr, Xr), 

PQPSr\QPXslQPx,.lSr:QPYr\Xr,Yr,Sr.Q-\'^\-^ 

I{Sr;Yd\Q) + I{X„Xr;Yd\Sr,Q) - I{Yr;Yr.\X,,Xr,Yd,Sr,Q)} 

> r(NNC-HD) C4il, (43) 



Appendix F 
Achievable rate with NNC 

The largest achievable rate for the memoryless relay channel 
is the combination of PDF and CF/NNC proposed in the 
seminal work of Cover and ElGamal ID. Here we use NNC 
to show: 

Proposition 12. The capacity of the G-HD-RC is lower 
bounded as 

^(HD-RC) > ^(NNC-HD) 

:= maxmin |4^^^'' + ^ Jijh.ij, 7»jAo,»j }; 



(jj)e[0:l]2 



where the maximization is over 



(i,i)e[0:l]2 



1, 



e [0, 1] : 

(»j)e[0:l]2 

Ps;, > : P-'^ ^ 1' 

(*j)e[0:l]2 



(»j)e[0:l]2 



(44a) 

(44b) 
(44c) 
(44d) 



where the different mutual information terms in ( |44| l ore 
defined next. 

Proof The NNC scheme in [4] Remark 18.5] adapted to 
the HD model gives the rate lower bound in (|43| at the top of 
this page, where the mutual information terms {IgAj,Iivi.ij}, 
€ [0 : 1]^ and in ( |44a| i are obtained as follows. 

We consider the following assignment on the inputs and on 
the auxiliary random variables for each G [0 : 1]^ 



= i, Sr = i] = ^ij such that ( |44b[ ) is satisfied 



Q=i,Sr-=3 



'-3,1 





r,i] 



such that \AAc\ and ( |44d[ ) are satisfied, 

Yr\x^,Y,^,Q=i,S,-=j = Yr + Zr^ij , 



A/^(0,(T? ) and independent of everything else, 



and in order to meet the constraint that Xg cannot depend on 
Sr conditioned on Q we must impose the constraint that in 
state Q — i, Sr — j the power of the source only depends on 



the index i. Then for each € [0 : 1]^ 

IiX,;Yr,Yd\Xr,Q = l,Sr=j) 

- log + + TT^) ) 

I{X„Xr;Yd\Q = l,Sr=j) + 

- I{Yr; Yr\Xs,Xr, Yd, Q=l, Sr = ]) 

= log(l + 5F,,, + /jF,,„)-log {^ + ^ - ^9 
I{Sr; Yd\Q) = ~Y '^iJ log(^ii) - log(7re) 
log- 



(45) 



(46) 



+ (700 +701) E 
+ (710 +711) E 

+ 

YdlQ^i ^ flit) := 



/„(r)IO = » 

log 10^1 

700 1 



(NNC) 







7oo + 7oi •^^'00 
701 1 



700 + 701 TTWoi 
710 1 



710 + 711 '^^'10 
7ii 1 



710 + 711 """ii 



exp(-|i|Vwoo) 
exp(-|t|Vwoi), t e 
exp{-\tf/vio) 
exp(-|t|Vvii), t € 



:= Ya.r[Yd\Q ^ i,Sr = j] = I + S Psa + I j Pr,ij- 
This proves the lower bound in (|44]) as a function of 

4,(z,j)e{0,l}2. 

In order to find the optimal af^, {i,j) G {0, 1}^ we reason 
as follows. /io,ij in ( |45] l is decreasing in afj while /g ^ in 



( [46| ) is increasing. At the optimal point these two rates are the 
same. Let 

CPs A 



1 + SPs. 



r ■.^I{Sr,Xr;Yd\Q), 
and rewrite the lower bound in (|44]) as 

r(^''^-™) = (7oo + 7oi)log(l + 5P,,o) 
+ (7io+7ii)log(l + 5'P,,i) 
- 700 log (1 + xq) - 710 log (1 + xi) 

+ min I700 log (I+xoCq) + 710 log (1+xiCi) , /'|. 
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The solution of 

min |7oolog(l + Xo) +7iolog(l + [ 
subject to 700 log (1 + xqCq) + 710 log (1 + xid) = /' 
can be found to be 

(1 - v)Ci 

with 77 < 1 such that 

700 log (1 + xoCo) + 710 log (1 + xiCi) = /'. 

■ 

Remark 5. For the special case of Q — Sr, that is, 
^ j(^Sr;Yd\Q) = I{Q;Yd\Q) = 0, the achievable 
rate in Proposition [72] reduces to 

^(NNC-HD) > l^^g ^ _ ^y^^^ 



(7,/3)e[0,l]" 

7/11 + (1 - 7)-^12|, 



/9:=log(l + 5P,,o)-log 1 + ^ 

V 0-0 

/lO :=log(l + 5P,u+/Pr,i), 

/ll ~ log (1 + SPs.Q + T-^Ps,Q 

\ 1 + '^0 

/12 :=log(l + 5P,,i). 
B + 1 



(1 + ^)^-^-1' 



A :-- 



IPr,l 






1 + SPs.l ' 


1 


P 

— —1 -Ts.l 

7 


1 
1 


- 7 



1 + 5P,,o ' 
Pr,l = 



1-7' 



(48a) 

(48b) 
(48c) 
(48d) 
(48e) 
(48f) 

(48g) 
(48h) 



where the optimal value for (Tq in ( |48f| l is obtained by 
equating the two expressions within the min in (|48a|l. 



Proposition 13. NNC with deterministic switch achieves the 
gDoF upper bound in (|3]l. 

Proof: With the achievable rate in Remark |5] (where here 
we explicitly write the optimization wrt ctq) we have ( |47] i at 
the top of next page, where we defined C5 and cg as 

C5 = ci := — ^ ^ ^ — --r^ > 1 since / > and as in (|27l 
log (1 + 5) 

log (1 + ^+5) 
C6 1 , ° , ^ > 1 since C > 0, 



log (1 + 5) 



and where 

7nnc := 7 TY^TT ^ [0,1J. 

(C5 - 1) + (cg - 1) 

By reasoning as for the PDF in Appendix |Dj it follows from 
the last rate bound that NNC also achieves the gDoF in (|3]l. 



Remark 6. For the special case of Q = fb, i.e., the time- 
sharing variable Q is a constant, the achievable rate in 



Proposition 12 reduces to 

^(NNC-HD) 



max mini I \Xs]Yr,Yd\SrTXrl 



I {Xs ; Xr ,Sr',Yii) — I ;Yr\Sr, X^ ,Xs,Yiij | 

c 

mm < 7 log I i + 5H — 

7e[o,i],T 



> max min<{7log ( 1 + 5+ ^ " _^ ) +(1-7) log (l + S") , 



1 + ^2 



/(5,;yd)+7log(l + 5)-7log 1 + 



+ (l-7)log 1 + 5 + 



1-7 



Note that with Q = $ the source always transmits with 
constant power, regardless of the state of the relay, while 
the relay sends only when in transmitting mode. Thus in this 
particular setting there is no coordination between the source 
and the relay. 



Appendix G 
Proof of Proposition[7] 

Proof: Consider the upper bound in (|5]l and the lower 
bound in (|8]l. Since the term I{Xs,Xr, 5,-; Yj) is the same in 
the upper and lower bound, the gap is given b}j^ 

GAP <IiXs;Yr, Yd\Xr, Sr) - I{U- Yr\Xr, Sr) 

-I{Xs;Yd\Xr,Sr,U) 

Next we consider two different choices for U : 
• For C < 5 we choose U — Xr and 

GAP < I{Xs;Yr,Yd\Xr, Sr) - I{Xs;Ya\Xr, Sr) 

= I{Xsi Yr\Xr, Sr, Yd) 

= V[Sr = 0]I{X,;VCX, + Zr\Xr, Sr = 0,\/5X, + Zd) 

+ V[Sr = l]I{Xs; Zr\Xr, Sr = 1, ^/SX^ + Zd) 

< V[Sr = 0] I{X^; VCXs + Zr\VSX, + Zd) 

< ¥[Sr - 0] log(l + C/il + S)) 
<log(l + 5/(l + 5)) 

< log(2). 

^ Let a lower bound be min_4{/;(y4)} and an upper bound be 
min^{/ii(^)}. With the definition 

^u,min := argmin{/„(^)}, A,min := arg min{/i (yl)}, 
A A 

we have fu {J^u min ) ^ fu{A.i min ) • This fact implies that the gap is upper 
bounded as 



GAP < min{/„(^)} - min{/,(^)} = /n(X.min) - /i(A.min) 
A A ■ 

< /n(^i,min) - /!(^!,min) < max{/„(^) - fi{A)}. 

A 
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^(NNC-HD) > 



+(1 - 7) log 1 + 



7e[0,l],<Ti^>0,/3G[0,l] 

(l-/?)5 



, , , PS\ , / 1 

mm < 7 log 1 H 7 log IH ^ 

7 / \ 



1-7 



/ 

1^ 



7 log 1 



C/3 



*^V(l-,)log(l+<lj^) 



(l + cT^h 7 

/3=7 

> max min{7log (1 + S") + (1 - 7) log(l + 5" + /) , 
7e[o,i],T^>o 

^ + ^ V (1 - 7) log (1 + - 7 log f 1 + A 



7 log 1 



1 + 



= max 
7e[oa].<T2>o 



log(l + 5)min{7+ (1 - 7)05,706 + (1 - 7)} - 7 log 1 + 



<t2>0 \ (C5 - 1) + (C6 - 1) 

t'-log(2) + log (1 + ^)^1+ 



log(l + A 



log 1 



c 



(1+<t2)(1 + S) 



(C5 - 1) + (C6 - 1) 



(47) 



• For C > S we choose U = X^Sr + ^^(l — 5*,^) and 

GAP <I{Xs;Yr,Yd\Xr.,Sr) 
- J(X«;yd|X,, 5,, XrSr + Xs(l - Sr)) 

= V[Sr = 0] (/(X,; y„ Ya\Xr, Sr = 0) 
-/(X,;y,|X,,^, =0)) 

+ V[Sr = 1] (/(X,; Yr, Yd\Xr, Sr = 1) 
-I{X,;Yd\Xr,Sr = l)) 

= V[Sr = 0] IiXs;Yd\Xr, Sr = 0, r,) 
+ ¥[Sr = 1] /(X,;y,|X,,,5, = l,Ya) 
= V[Sr = 0]I{Xs; y/SXs+Zd\Xr, =0,\/CX, + Z,,) 

+ F[Sr = 1] /(Xs; ^rl^r, Sr = 1, + Zrf) 

< PiS*^ = 0] liXs; VSXs + ZdlVCXs + Zr) 

< ¥[Sr = 0] log(l + S/il + C)) 
<log(l + C/(l + C)) 

< log(2). 

This concludes the proof. I 



Appendix H 
Proof of Proposition[8] 



Consider the upper bound in (|6]l and the lower bound in (|9|. 
Recall that Ii — I5 I2 — h I3 > I7 I4 ^ h and therefore 



Proof: 



GAP < max 



{^(7) + 7^1 + (1 - l)l2 - ih - (1 - 7)^6, 

7/3 + (1 - 7)/4 - ih - (1 - 7)4} 
= max {^(7), 7(73-/7)} 

< inax< log(2),log ^, ^ > 

1^ \l + max{C,S}PsfiJ i 

< max <^ log 2), log —— \ 

I \ 1 + maxjC, 5}-Ps^o / J 

< max I log(2), log(2)| log(2) = 1 bit. 

This concludes the proof ■ 

Appendix I 
Proof of Proposition[9] 

Proof: Consider the upper bound in (|7]i and the lower 
bound in ( fTTj i. We distinguish two cases: 

. Case 1: 5 > C. hi this case r^^DA-HD) ^ j^g^-j^ ^ gy 
The gap is 

GAP < r^^^"™) - r^^^^-™} 

(61 - 1)(62 - 1) 



<21og(2) + log(l + 5) 
<21og(2) + log(l 
= 21og(2) + log (1 



ih - 1) + {b2 - 1) 
<21og(2) + log(l + 5) (62 -1) 
C 



< 21og(2) + log i^l 

< 31og(2) = 3 bits. 



1 + S 

S 
l + S 



• Case 2: S < C. First, by noticing that 

log (1 + (%/7 + y:S)2) < log (1 + / + 5) + log(2), we 
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further upper bound the expression in (|7]l as 

^(CS-HD) <21og(2)+log(l + 5) 

log(l + T^) + log(2)) log (l + ^ 
log ( 1 + ^) + log(2) + log (l + 



Next we further lower bound r'-^^^ hd) jj^jj 

^(LDA-HD) >iQg(^^^) 

l°g(l + TT5) (log(l + Tf5)-log(2) 



log 



(l + ^)+log(l + ^) 



Hence, with x = log ( 1 + j ,y = log ( 1 + j , 



we have 



GAP < r'^'^S-HD) _ ^(LDA-HD) 

{x + l)y _ x{y - 1) 
X + 1 + y X + y 



< 2 



= 2 



x'^ + y'^ + xy + X 



10, 



a;^ + j/^ + 2a;j/ + a; + ?; 
< 3 bits. 

This concludes the proof. 

Appendix J 
Proof of Proposition[To1 

Proof: With NNC we have 

GAP < max |-H(7) + 7/1 + (1 - 7)12 - 7^ - (1 - 7) 

7/3 + (1 - - "/hi - (1 - 7)-^i2| 
= max 1^(7) + 7log (1 + SPs^) + 7 log (^1 + 

+ (1 - 7) log (1 + {^/SI\i + VlB^if) 

- 7 log (1 + SP.^o) - (1 - 7) log (1 + 5P,a + IPr,i) , 
7 log (1 + (C + 5)F,,o) + (1 - 7) log (1 + SPsa) 

- 7 log (^l + 5P.,o + ^) -(I-7) log(l + 5P.a) } 
< max (^(7) + (1 - 7) log(2) + 7log f 1 + 

7 log 1 + 



1 + SPsfi + Y^^CPsfl 



< max |-H(7) + (1 - 7) log(2) + 7 log + ■ 

7log(l + CTo)} 

< 1.6081 bits. 

where for CTq we chose the value 

^,^^^pHM±(1^7)log(2) 



by equating the two arguments of the max (this is so because 
"^(7) + (1 - 7) log(2) + 7log (^1 + is decreasing in CTq, 



while log (1 + ctq) is increasing in ctq). Numerically one can 
find that with the chosen ctq the maximum over 7 G [0,1] is 
1.6081 for 7 = 0.3855. 

Note that by choosing 0-0 = 1 the gap would be upper 
bounded by 2 bits. ■ 

Appendix K 
Proof of Theorem|2] 

Proof: The capacity of the HD channel in (|29| is upper 
bounded by the capacity of the FD version of the same 
channel, which is given by 



= max mm\l{Xs,Xr;Yd),I{X,;Yr,Yd\Xr)} 
= max WLUiiH{Yd),H{Yr,Yd\Xr)\ 

PXs,X,. < J 

I max{/3sd, /3rd}, max{/3sd, /3sr}| 



(a) 

= mm • 



= /33d + min{[A.d - ^sd] + , [/3sr - /3.d] + }. 



(49) 



where the equality in (a) is in general an upper bound but 
for this channel model is achieved with equality with i.i.d. 
Bemulli(l/2) input bits. 

For the capacity of the HD channel we distinguish two 
cases: 

. Regime 1: /3rd < /3sd or (3^^ < /^^d- 

In this regime, C^hd) < c'(fd) ^ gjjj^g ^^le rate 
can be achieved by silencing the relay and 

using i.i.d. Bernulli(l/2) input bits for the source, we 

conclude that (7*^^^' — (3sd in this regime. 
. Regime 2: ^,.d > Psd and (5^^ > (3sd- 

We start by writing = [Yd^u,Ydd], where 

- Yd^i contains the lower /3sd bits of Yd- These bits are 
a combination of the bits of Xg and the lower bits of 
Xr- The lower bits of X^ are indicated as X,. With 



reference to Fig 2(b) Yd^i corresponds to the portion 



of Yd containing the "orange bits" labeled 61 [2]. 
- Yd^u contains the upper f3^d — P^d bits of Yd- These 
bits only depend on the upper bits of Xr- The upper 
bits of Xr are indicated as X^ „ ■ With reference 



to Fig 2(b) Yd^u corresponds to the portion of Yd 



containing the "green bits" labeled a. 
We have 

H{Yd)=H{Yd,u,Ydd) 

< H{Yd.u) + H{Yd,i) 

<H{Yd,a) +Psd, 

since Yd,i contains jS^d bits and where H{Yd,u) is com- 
puted from the distribution 

ViYd.u ^y]= = 0]P[rd,„ = y\Sr = 0] 

= 7(5[?/] + (1 - 7)P[^r.« =v\Sr = 1] 

for e [0 : L - 1], L := 2P-'^-f^-^ > 1, where 5[y] = 1 
if y = and zero otherwise, and where 7 := P[S'.r = 0]. 
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max{R} < I{Xi,[XA,SA];YA^,YK\[XA^,SA^], Si = 1,Sk = 0) 

< H{Sa\Sa^, 5i = 1,Sk - 0) + I{Xi,Xa;Ya^,Yk\Xa^,S[2..k-i], - 1,^a' = 0) 



(b) 



< HiSA)+ J2 A, log|l|^e|+i+H^.,K 



{l}uA,stiA,s\ 



s=0 



s=0 



(c) ^ 

<|^|log(2)+ ^ A,Rank[H^,,]log(max{l,Trace[K{i}u^,.]}) 

s=0 

+ A. log|l|^=|+i +h^,,h5_,| 

s=0 

< 1^1 log(2)+min(l + 1^1,1 + 1^^1) log max i ^ A„ ^ A, Trace[K{i}uA,] 



s=0 



s=0 



J2 A. log|l|^=|+i +H^,,H5_, 



(o) 



< 1^1 log(2)+min(l + 1^1,1 + 1^^1) log (1 + 1^1)+ ^ A, log |l|^e|+i + H^,,H5_, 



(50) 



Let ¥[Xr,u = y\Sr = I] = Py e [0,1] : ^^^p^, = 1. 
Then 

i?(l^d.«) 

= f([7 + (1 - 7)po, (1 - 7)pi, . . . , (1 - 7K-1]) 



< ff([7+(l-7)po,(l-7) 



1 -Po 



.(1-7) 



1 -Po- 
L- 1- 



(l-^?)log 



1 



L- 1 



The upper bound on the entropy H{Yd «) is maximized 
by 



= 1 - max{l/L, 7} <+=+ pp = 



[I/L-7]- 



Therefore we have 

(7(HD) 

im{l{Xs,Xr,Sr;Yd),I{X,;Yr,Yd\Xr,Sr)} 

ax mm{H{Yd),H{Yr,Yd\Xr,Sr)} 

< /3sd + max min | (1 - 6**) log + 0* log ^— ^ 

-ys 0,1 L 1 — t^* 6^* 



= max mm • 

PXs,Xr,Sr 



max 



7[/3sr - A. 



H{Yr, Y^lXr, Sr) 

= ¥[Sr=Q]H{Yr,Yd\Xr,Sr=0) 
+ ¥[Sr = l]H{Yr,Yd\Xr,Sr = 1) 

< 7max{/3sr,/3sd} + (1 - 7)/3sr- 



In order to show the achievability of pOf consider the 
following input: the state Sr is Bernulli(l — 7) indepen- 
dent of everything else and Xg and X,, are independent. 
The source uses i.i.d. Bernulli(l/2) bits. The relay uses 
i.i.d. Bemulli(O) bits for Xr.i and P[XrM ^ v] ^ Pq if 
J/ = and P[Xr;u = y] = (1 - Po)/iL - 1) otherwise, 
i.e., the components of Xr.u are not independent. Notice 
that the distribution of Xr.u in state Sr = is irrelevant 
because its contribution at the destination is zero anyway, 
so we can assume that the input distribution for Xr is 
independent of the state Sr- 



Appendix L 
Proof of Theorem|4] 

Proof: 

a) Upper bound: The cut-set upper bound on the capac- 
ity of the HD Gaussian relay network gives, for each A, ( jSO] ) 
at the top of this page where the inequalities are due to the 
following facts: 

• Inequality (a): chain rule of the mutual information; 

• Inequality (b): by considering that the discrete random 
variable S[2-k-i] has at most 2^^^ masses and by letting 
As P[5[2:K-i] = .5] e [0, 1] for s e [0 : 2^-^ - 1] 
such that X]s=o As = 1. 

Here we use the convention that "5[2:A'_i] ~ s" means 
that the j-th entry of S[2:K-i] is equal to the j-th digit in 
the binary expansion of the number s. For example: with 
= 5 and s = 4 = 1 • 2^ + • 2^ + • 2°, the notation 
"S[2:K-i] = s means ^3 = 1,^3 - 0,^4 = 0". K^,^ 
represents the covariance matrix of Xa conditioned on 
[S[2:K-i] ~ Si — 1,Sk — 0] and liA,s is the matrix 
obtained from (I — S)HS by retairung the rows indexed 
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iR&i,{R} > I{X^,XA-,YA^,YK\XA^,XK,Sy2-,K-i\.Si = 1,Sk = 0) - I{Yx:YA\Xy^-.K\,YA^,^^ = 1, ^if = 0) 



> J2 A, log I|^e|+i + — ^ H^^.H^, -|^|log(l , 2 

> J2 A, log|l|^e|+i+H^,,Hi,| + J2 A, Rank[H^^,] log ( mini 1,-^1) -1^1 log + ^ 

2^-2-1 

> ^ A, log |l|^e|+i + H^,,Hi,| -min(l + 1^1,1 + 1^=1) log(l + a2)- 1^1 log [l + ^ 

2^~^— 1 

V' A, log|V|+i+H^,,H5_,|-inin(l + 2|^|,X-l)log(2) 



(51) 



s=0 



by {K} U and the columns indexed by {1} U ^ for 

diag[S] = S[i.K] such that [S[2-k-i] = s, Si ^ 1,Sk = 
0]. 

Because of the power constraint we must have 
E'1V"'aJK[i:K]J <1, ke[l:K]. 

• Inequality (c): by exploiting the following relation: ^ 
K ^ Amax(K)I ^ Tracc[K]I. Moreover for a 7^ and 
by using the eigen-decomposition K = UAU^ G C" 
with [A] 11 > [A] 22 ■ • ■ [A]„„ the following holds 

ii+H Ki=ii+H A|= n 

j" = l,...,Rank[K] 

<max{l,|a|}^^"'^[Kl [] {I + Xj) 

j = l,...,Rank[K] 

= max{l,|a|}^^"'^[^l |I + K|. 

• Inequality (d): since the rank of a matrix is at most the 
minimum between the number of rows and columns. 

• Inequality (e): since the entropy of a discrete random 
variable is at most the log of the cardinality of its support 
and because of the input power constraints. 

b) NNC lower bound: A lower bound to the capacity of 
the memoryless HD-MRC is found by adapting the NNC for 
the general memoryless MRC ifTTI to the HD case. In all states 
we consider i.i.d. Af {0,1) inputs with Q = S^i-x] and with 
Yk i^fc + ^fc for Zfe ~ A/'(0, CT^) independent of everything 
else. With this, the NNC lower bound gives, for each A, ( |5T| ) 
at the top of this page. 

c) Constant Gap: The gap between cut-set upper bound 
and the NNC lower bound is 



GAP < max 

\A\e[Vi.K-2 



[\A\ log(2) +min(l +2|^| , if - 1) log(2) 

+ min(l + \All + \A''\) log (1 + \A\) } 
max { min(l + \Al 1 + log (1 + \A\) 

min(l + ?,\A\, \A\+K- 1) log(2)} 



l-4|e[o 



Appendix M 
Proof of Proposition[TT] 



Proof: The proof of Proposition 1 1 follows directly from 



the proof of Theorem |4j by taking into consideration the 
following facts: 

• In the diamond network the channel matrix H has rank 
2. Thus Rank[H_4 s] — 2 both in the cut-set upper bound 
and in the NNC lower bound; 

• Instead of the upper bound H{Sa) < \A\ log(2) we 
proceed as follows 



'2^-^-1 

s=0 



{A,>0} 



where l{;^_,>o} is the indicator function defined as 



L{A,>0} 



if > 
otherwise 



Let A* be the optimal values from the cut-set upper bound. 
Then the gap is given by ( |52j i at the top of next page. The 
optimal value of the quantization noise is <t^ = K/2 — 1 and 
it has been found by equating to zero the first derivative of 
/(ct2) = 21og(l + cr2) + (K -2)\og{l + Now to be 

sure that this cr^ is the one that minimizes /(cr^) we need to 
compute the second derivative that is given by 

-2cr^ + {K - 2){2a^ + 1) 



The denominator of f icP') is always a positive number Now 
we need to evaluate the numerator of f"{(J^) in cr^ = K/2 — 1 
and verify that it is always positive. We obtain 

2 



2a^ + {K- 2) (2(7^ + 1) 



K-2 



-{K~2){K-1) 



= K 



K 



- 1 > V if > 3. 



This shows that for a diamond network the optimal quanti- 
zation noise is cr^ — K/2 — 1 and that the gap has two 
components: one depends on how many non-zero AJ,s G 
[0 : 2^~^ — 1], are needed to attain the cut-set upper bound 
(which in general is upper bounded by 2^~^), and the other 
is logarithmic in the number of nodes in the network. 
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GAP<min^^^m;«_^Jlog 1^ ^ l{A*>o}j + 21og(l + |^|) + 21og(l + a^) + |^| log(l + l/a^)} 
= l°g| E 1{A*>0} I +21og(if-l)+min{21og(l+a2) + (i^-2)log(l + l/a2)} 

^ l{Aj>o} J +21og(if-l) + 21og(^l + :^^) +(if-2)log(^l + -^) 



<log ^ l{A*>o} +21og(if(if-l)) +21og(e/2) 



5=0 



<{K-2) log(2) + 41og(i^) + 21og(e/2). (52) 



Recently, it has been shown numerically for diamond net- 
works with at most seven relays |17|, i.e., K < 7, that 
X]s=o l{A*>o} 1^ if — 1. In the same work it has been 
conjectured that the same holds for any number of relays. 
If the conjectured is true then our gap would be at most 
51og(i^) + 21og(e/2). 
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Fig. 1: The Half-Duplex Gaussian Relay Channel. 
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(a) HD Phase I {Sr = 0). 



(b) HD Phase II (Sr = 1). 
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Fig. 2: The gDoF optimal strategy for the linear deterministic approximation of the Gaussian noise channel at high SNR. HD 

in Figs. |2(a)| and |2(b)| and FD in Fig. |2(c) 




3: Comparison of the capacities of the LDA for both HD and FD. Also, upper and lower bounds for the capacity of the 

HD channel. 
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Fig. 4: Comparison of the rates of the G-RC without a direct link for 5 = 0, C = 15, / = 3. 
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5: Comparison of the rates of the G-RC for 5 = 30 dB, C = 37.63 dB, I = 34.77 dB. 
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6: Comparison of the rates of the G-RC for S" = 30 dB, C = 32.55 dB, I = 40.41 dB. 
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Fig. 9: Numerical evaluation of the maximum gap varying the SNR for f3sd = 1 and {(3^^, Psr) G [lj2.4] with deterministic 

switch (red curve) and random switch (blue curve). 
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Fig. 10: Gap in ( [37] i (blue curve) and limit behaviour in ([38| (red curve) as a function of K with — 1. 
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a > 0, c e [0, 1], fo>0 

•^^=4 = inm{l + a, max{l + c, 1 + b}} > 1 

<^=i.bssueiay = "^infl, m^{l + a, (1 + c)/2}} = 1 

Fig. 11: Example of a full-duplex 2-relay network with gDoF strictly larger than the gDoF obtained by using the best relay 
only. The numerical value on a link represents the SNR exponent on the corresponding link. 



